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The maximally entangled state can be in a mixed state as well as the well-known pure state. Taking the 
negativity as a measure of entanglement, we study the entanglement dynamics of bipartite, mixed maximally 
entangled states (MMESs) in multipartite cavity-reservoir systems. It is found that the MMES can exhibit the 
phenomenon of entanglement sudden death, which is quite different from the asymptotic decay of the pure- 
Bell-state case. We also find that maximal entanglement cannot guarantee maximal nonlocality and the MMES 
does not correspond to the state with maximal measurement-induced nonlocality (MIN). In fact, the value and 
dynamic behavior of the MIN for the MMESs are dependent on the mixed state probability. In addition, we in¬ 
vestigate the distributions of negativity and the MIN in a multipartite system, where the two types of correlations 
have different monogamous properties. 
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I. INTRODUCTION 

The maximally entangled state plays an important role in 
quantum information processing (ll-l^, including quantum 
teleportation |3l, quantum cryptographic protocols i] and 
quantum dense coding |@]. In bipartite d systems, Caval¬ 
canti et al proved that all maximally entangled states are pure 
states a. For example, in the simplest 2 0 2 systems, the 
pure maximally entangled state is the Bell state which can be 
written as 

|^) = (|00) + |11))/V2 (1) 

up to local unitary transformations. Recently, it was fur¬ 
ther shown that there exist mixed maximally entangled states 
(MMESs) in bipartite d ® d' systems with d' > 2(i, which 
can be used as a resource for faithful teleportation The 

MMES in 2 (g) 4 systems has the form |H] 

P = P|V’1)('*/'1| + (1-P)IV’2)('*/'2|, (2) 

where the mixed state probability p lies in the range (0,1), and 
the two pure state components are = (|00) -I- \ll))/y/2 
and IV' 2 ) = (|02) + |13))/y2. 

The dynamic behavior of entanglement is a fundamental 
property of quantum systems. This is because unavoidable in¬ 
teractions with the environment may lead the entanglement of 
quantum systems to be degraded and, in certain cases, to dis¬ 
appear in a hnite tim^i.e., the so-called entanglement sud¬ 
den death, ESD) IUOI - It^ . Lopez et al analyzed the dynamic 
behavior of entangled cavity photons being affected by two 
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dissipative reservoirs ifl^ and found that the entanglement of 
cavity photons initially in the two-qubit Bell state decays in 
an asymptotic manner. However, for the newly introduced 
MMES, its entanglement dynamic property is still an open 
problem, especially for a real quantum system. Since the 
MMES is a perfect physical resource in quantum information 
processing 0,1^, it is desirable to investigate its dynamical 
property in a quantum system. This is because, once the en¬ 
tanglement evolution experiences the ESD, we are no longer 
able to concentrate the entanglement of MMES, which results 
in some entanglement-based quantum communication proto¬ 
cols losing their efficacy. In this sense, the study of the en¬ 
tanglement dynamic property of the MMES can provide not 
only useful knowledge for practical quantum operations but 
also the necessary information to cope with the decay of en¬ 
tanglement. 

Nonlocality is also a kind of resource in quantum informa¬ 
tion processing lll^ a nd has a close relationship with quan¬ 
tum entanglement iflal . The measurement-induced nonlocal¬ 
ity (MIN) ifl^ is the maximum global effect caused by lo¬ 
cally invariant measurement, which is different from the con¬ 
ventionally mentioned quantum nonlocality related to the vi¬ 
olation of Bell’s inequalities lEol 1211] . Moreover, the MIN 
can quantify the nonlocal resource in quantum communica¬ 
tion protocols involving local measurement and a comparison 
between the pre- and postmeasurement states S. Luo and 
Eu proved that, for the pure Bell state, the MIN achieves the 
maximal value m. But it is not clear whether or not the 
MMES also has the maximal nonlocality. In particular, can 
the maximal entanglement guarantee the maximal nonlocal- 
ity? Eurthermore, in order to obtain a deep understanding of 
the dynamic properties of MMESs, it is helpful to analyze the 
entanglement and nonlocality distributions in an enlarged sys¬ 
tem including its environment. 

In this paper, as quantihed by entanglement negativity 
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and the MIN Cl, we study the dynamic properties of the 
MMES in the dissipative procedure of multipartite cavity- 
reservoir systems. It is found that the MMES can disentangle 
in a finite time, which is quite different from the asymptot¬ 
ical decay of a pure Bell state. We also find that the MIN 
of the MMES is not maximal, which means that the maximal 
entanglement cannot guarantee the maximal nonlocality. In 
addition, the evolution of the MIN is dependent on the mixed 
probability of the MMES. Einally, we investigate the distribu¬ 
tions of the negativity and the MIN in the multipartite system, 
where the squared negativity is monogamous but the MIN is 
not monogamous. 


II. DYNAMIC PROPERTIES OF ENTANGLEMENT AND 
NONLOCALITY FOR THE MMES 


We first recall the definition of the MMES before analyzing 
its dynamic properties. Ind^ d' systems, a mixed state is an 
MMES if and only if it has the form |H,|^ 

K 

P~ ^ ^ Pm IV^m) ('^m I, (3) 

m—1 

where the mixed state probabilities satisfy X]m=i P™ = 1 
with K < floor(d'/c?), and the pure state component is 


IV'm) 


1 

7 ^ 


d-l 

|i) 0 |i -f (to - l)d), 

i=o 


(4) 


which is the maximally entangled pure state. In the following, 
we will study the dynamic properties of MMESs in bipartite 
2 (g) 4 systems. 

We consider a practical dynamic system of two cavities in¬ 
teracting with two independent reservoirs. The initial state of 
the four-partite system is 


PciC2r-ir2 (0)=pciC2(0)®|00)(00| 

rir 2 (5) 

where the two reservoirs are in the vacuum state, and the two 
cavities are in an MMES, 

PciC 2 ( 0 ) =p|V'l)(V'l| + (1 - 4 ')IV’ 2 )(V’ 2 | ( 6 ) 

with = (|00)-f|ll))/V2and IV’ 2 ) = (|02)-f |13))/V2. It 
should be noted that, although PciC 2 is written as a probability 
mix of IV'i) and \ip 2 ), its pure-state component has the generic 
form ■yyql'tpi) + e'^’^yyi — q\il> 2 ), with the parameters q G [ 0 , 1 ] 
and (j) G [0, 27 r]. The interaction of a single cavity and an N- 
mode reservoir is described by the Hamiltonian lllal23l4^ 

N N 

H = hua'd -t- h'^^uJkb\bk + h'^^gk{db\ + bkd))- (7) 

fc=l k=l 


At later times, in the limit N ^ 00 for the reservoirs with a 


flat spectrum ifl^ , the state is given by 

P t t 

f’ciJ-iC2r2 (i-) ~ '^{\4’o)ciri C2r2 \4*l)ciri\4^l)c2r2) 

x((</> 0 | ciri (M C2r2 + (<(’l|ciri('('l| C2r2)] 

+ 2 c^rx\4'2) C2r2 + I'/'l) ciri j (/>3)c2r2 ) 

x((</> 0 | ciri ((/> 2 |c 2 r 2 + (^l|ciri((/>3|c2r2 )], ( 8 ) 

where the components can be written as 

= |oo), 

I7i) = C| 10 )+x| 01 ), 

14) = e'| 20 ) +72^x111)+x"| 02 ), 

14) = 4|30) + 434x|21) + 43ex"|12) + x'|03),(9) 

in which the amplitudes are and x(i) = (1 ~ 

e“'^‘)i/ 2 , with the parameter k being the dissipative constant 

M- 

In this section, we focus on the dynamic properties of two 
cavities which are initially in an MMES. As the cavities and 
reservoirs interact, the state of two cavities is PciC 2 it) = 
Trrir -2 [PciriC 2 r 2 (^)], which has the matrix form 

f ctii 0 0 0 0 uiQ 0 0 ^ 

0 U 22 0 0 0 0 U 27 0 

0 0 033 0 0 0 0 038 

, , _ 0 0 0 044 0 0 0 0 

PciC 2 (t)- 0 0 0 0 055 0 0 0 

qqi 0 0 0 0 qqq 0 0 

0 a72 0 0 0 0 0)77 0 

\ 0 0 083 0 0 0 0 088 / 

( 10 ) 

with the basis in the order 

{|00),|01),|02),|03),|10),|11),|12),|13)} and the ma¬ 
trix elements 

an = {P + X^ + X^ - PX®)/ 2 , 

022 = 4x^[2 - p -f 3(1 - p)x'‘]/2, 

033 = (1 -p)4(l + 3x'‘)/2, 

044 = (1 -p)C®xV2, 

055 = 4x7p + X^-PX^)/2, 

066 = 4b+ 3(1 -p)x'‘]/2, 

077 = 3(1 -p)^®xV2, 

088 = (1 -p)4/2, 

016 = 061 = 4b + 73(1 -p)x'‘]/2, 

027 = 072 = \/72(1 - 44x^, 

038 = 083 = (1 — p)C®/2. (11) 


In order to characterize the dynamic entanglement proper¬ 
ties of two cavities, we need to choose a suitable measure of 
entanglement. Here, we use the negativity iIt^I to quantify the 
entanglement of two cavities due to its computability for any 
state of an arbitrary bipartite system. Eor the quantum state 
PciC 2 it), its negativity is 


iy=ic 2 (t) — 


llpl4(i)ll-i_ElilAd-i 


2 


2 


( 12 ) 
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FIG. 1: (Color online) Evolution of entanglement of the MMES in 
two cavities, where the negativity is shown as a function of the time 
hit and the probability p. The red line indicates the entanglement 
sudden death of the two cavities. 


where 11 • 11 is the trace norm and equal to the sum of the moduli 

T 

of the eigenvalues for the Hermitian matrix /9ciC2 {t), which is 
the partial transpose with respect to the subsystem ci Q. 
After some derivation, we can obtain the eigenvalues 

Ai = (i-p)eV2, 

>^2 = {p + x^ + - PX®)/2, 

As = ?"{[! +(6-6p)x"]-^}/4, 

A4 = ^"‘{[1 + (6 - 6 p ) x ‘^] + \/ A }/ 4 , 

A5 = [2(l-p)CV-\^]/2, 

A6 = [2(l-p)eV + v^]/2, 

Xr = e{x^[l + 2il-p)x^]-VC}/2, 

A8 = C"{x"[1 + 2(1-p)x^] + ^}/2, (13) 

where the parameters are A = (1 — 2p)^ + 24(1 — 

B = {1- + x‘^)> and C = + (1 - p)[l + (2v/3 - 

l)p]x'‘ + 5(1 - pfx^ + (1 - Pfx^"^- 

In Fig. 1, we plot the negativity NciC 2 {'t) as a function of 
the time nt and the probability p. As seen from Fig. 1, when 
Kt = 0, the quantum state PciC 2 is the MMES, and its nega¬ 
tivity has the maximal value of 0.5 regardless of the choice of 
probability p. For a given value of the probability, p, the nega¬ 
tivity decreases as the time, Kt, increases. It should be pointed 
out that as time increases, the entanglement of the MMES de¬ 
cays through sudden death rather than asymptotically like the 
two-qubit Bell state. The red line in Eig. 1 indicates the time of 


the ESD for the negativity satisfies the equation 

1_v/3 + 3x4-3x« + v^ 

V = -^- (14) 

l-2y3 + x-^ + 5x4-3x8 

where the parameter D = 3 — 2'/3 + 4(2 — '/3)x^ + X®> with 
X = \/l — e~'^* (the derivation of the ESD line is presented in 
Appendix A). When the probability p = 0, the initial state of 
the two cavities is \'ijj 2 ) = (|02) -F |13))/-\/2, which is a two- 
qubit pure maximal entangled state with qubit ci involving 
states 0 and 1 and qubit C 2 involving states 2 and 3. As the 
system evolves, the quantum state of two cavities becomes a 
2 0 4 system and has the matrix form shown in Eq.(lO) with 
the parameter p = 0. Its entanglement evolution shows the 
ESD phenomenon, and the negativity becomes zero at the time 
Kt = ln[(3 + •\/3)/2] ~ 0.8612. When the probability p S 
(0,1), the initial state PciC2(0) is the MMES. The ESD time 
of the two cavities is determined by Eq. (14) and increases as 
a function of the probability p. In the p = 1 case, the initial 
state is the two-qubit Bell state |'0i) = (|00) + \ll))/\/2, 
and its entanglement disappears at the time Kf —^ 00 , which 
coincides with the result for asymptotic decay presented by 
Lopez et al. S. 

Here, we have shown that, unlike the asymptotic entan¬ 
glement decay of the Bell state, the MMES of two cavities 
experiences the ESD in the dissipative procedure of cavity- 
reservoir systems. It is argued that the high-dimensional com¬ 
ponent IV' 2 ) = |02) + |13))/-\/2 plays the dominant role. Al¬ 
though the initial state 1 - 02 ) is a logic two-qubit state, it will 
evolve to a 2 (g) 4 system along with the cavity-reservoir in¬ 
teraction, which results in the ESD phenomenon of two cav¬ 
ities. In the evolution of two cavities, the ESD time is post¬ 
poned when the mixed probability of the component 101 ) (the 
Bell state) increases. In the case of p = 1 for the MMES, 
the component |02) disappears, and there is no ESD phe¬ 
nomenon, which is just the evolution of the pure Bell state 
|0i) = (|OO) + |ll))/V2. 

In addition to quantum entanglement, nonlocality is also 
a useful resource in quantum secure communication. It is 
worthwhile to further investigate whether maximally entan¬ 
gled states like the MMES also result in maximal nonlocality 
and how the nonlocality of the MMES evolves with time. The 
MIN {H is a computable nonlocality measure, which is the 
maximum global effect caused by locally invariant measure¬ 
ment. The MIN is defined as iS 

MIN(p^b) = maxnA||p- n^(p^B)|p, (15) 

where the max runs over all the von Neumann measurements 
H"^ = {life } which do not disturb the reduced density matrix 
PA (i-e., X)fe II^PaII^ = Pa), and the Hilbert-Schmidt norm 
is ||Ar|p = trX^X. The state of the two cavities PciC 2 {t) in 
Eq. (10) can be rewritten in a generalized Bloch form, 

15 3 15 

® E + E E ® ’ (16) 

j=l i=l j=l 
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where I 2 and I 4 are identity matrices of the subsystems and 
Xi = <Jil \/2 and Yj = ® crj ")/2 are operator bases with 

j'jj" G {0,1, 2, 3} except for the case j' = j" = 0 (here, 
(To = I 2 and {cti, (T 2 , (T 3 } are the Pauli matrices). In the Bloch 
expression, Eq. (16), the matrix elements are 

Xi — tr(pciC 2 ^t ^ -^ 4 / 2 ), 

Vj = tr(pciC2-^2/\/2(g)y,), 

Tij = tr{pc^c2Xi'S)Yj). (17) 

Luo and Fu derived an analytical formula for the MIN in an 
arbitrary 2 0 (i system El 

trTT^ ---^x*TT*x if x ^ 0 

II X IP 

trrT‘-A 3 ifx = 0(18) 

where A 3 is the minimum eigenvalue of the 3x3 matrix TT* 
with T = {Tij), and x = (a: 1 , 2 : 2 , 2 : 3 )* is the local Bloch 
vector with the norm ||x|p = (here, t represents the 

transposition). After some derivation, we can obtain the ex¬ 
pression for the MIN of two dissipative cavities, which can be 
written as 

MIN,,e,(f) = y^[F-2p{F-G)+p\l-2G + F)], (19) 

where the two parameters are F = + 3%® and 

G = '/Sx^ with ^ and % = \/\ — The details 

of the calculation and the continuity analysis of the MIN are 
presented in Appendix B. 

In Fig.2, we plot the MIN as a function of the time nt and 
the probability p. When nt = 0, PciC 2 ( 0 ) is the MMES, and 
its nonlocality is 

MIN,,,,(0)= (20) 

which is symmetric to the probability p = 1/2. As shown in 
Fig. 2, the MIN has the maximum value of 0.5 for the cases 
of p = 0 and p = 1 , which correspond to the pure maxi¬ 
mally entangled states 1 ^ 2 ) and lipi)- When p G (0,1), the 
MIN has less than the maximum value, and reaches the min¬ 
imum value of 0.25 at p = 0.5. Therefore, we can reach 
the conclusion that the nonlocality of the MMES is not maxi¬ 
mal, although its entanglement is maximal for any value of the 
probability p. According to Eq. (20), we find that the MIN is 
directly proportional to the purity Tr(p^) of the MMES, i.e., 
MINciC 2 ( 0 ) = iTi'[PciC 2 (0)]- When the mixed-state probabil¬ 
ity of the MMES changes from 0 to 1/2, its purity decreases, 
which results in the MIN changing from the maximum of 0.5 
to the minimum of 0.25. When the probability p changes from 
1/2 to 1, the purity of the MMES also increases, and the MIN 
changes from the minimum of 0.25 to the maximum of 0.5. 
As two cavities evolve, the MINciC 2 (i) decays in an asymp¬ 
totic manner and disappears in the limit Kt —^ 00 . This is 
different from the sudden-death evolution of the negativity of 
two cavities, since the nonlocality described by the MIN con¬ 
tains both quantum and classical correlations. The inset of 
Fig. 2 shows the difference MINp=i — MINp^g^ which indi¬ 
cates that the nonlocality is no longer symmetric as the system 
evolves. 



FIG. 2: (Color online) Evolution of the MIN of the MMES in two 
cavities, where the nonlocality is shown as a function of the time Kt 
and the probability p. The inset is the difference MINp^i — MINp=o 
as a function of Kt. 


III. ENTANGLEMENT AND NONLOCALITY 
DISTRIBUTIONS OE THE MMES IN MULTIPARTITE 
DYNAMICS 

Entanglement monogamy is an important property of mul¬ 
tipartite systems and means that quantum entan glem ent can¬ 
not be freely shared among many parties 1271 - 1^ . It has 
been proved that the squared negativity obeys the monogamy 
inequality in pure states of qubit systems, — 

^AiA 2 — ■ ■ • — ^AiAn — ^ However, for mixed states 
or multilevel pure-state systems, whether or not a similar 
monogamy relation holds is still an open problem. Recently, 
a numerical analysis was carried out for tripartite multilevel 
pure states 1^ . which supported the monogamy relations for 
squared negativity. However, it is still unknown whether or 
not the monogamous relation holds for the four-partite case, 
especially in a real quantum system with dissipative reser¬ 
voirs. With this in mind, we next analyze the negativity distri¬ 
bution of the MMES in the four-partite 2 ® 2 0 4 0 4 cavity- 
reservoir systems. On the one hand, this analysis can verify 
the monogamy inequality for the squared negativity, and, on 
the other hand, it can provide a deep understanding of the dy¬ 
namics of the MMES. 

The residual entanglement in monogamy inequalities can 
be used as a multipartite entanglement measure or indicator 
to characterize the structure of multipartite entanglement iIt^ 
[ 3 ^ . For composite cavity-reservoir systems, we analyze the 
entanglement distribution in the partition Ciri |c 2 r 2 and eval- 
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FIG. 3: (Color online) The distribution of negativity in the dissipa¬ 
tion of multipartite 2 (g) 2 (g) 4 (g) 4 cavity-reservoir systems, where 
the non-negative values of |c 2 r 2 indicate multipartite entangle¬ 
ment. 


uate the multipartite entanglement indicator 


initial state becomes the two-qubit Bell state. For the sub¬ 
systems Cir 2 and riC 2 , the negativities have two peak values 
when the probability is p = 0 and p = 0.5 [see Figs. 3(a) and 
3(b), where we multiply iVcjr 2 NriC 2 by a factor of 2 for 
clarity]. As the probability p increases, the number of peaks in 
the negativity changes from two to one, as shown in Figs. 3(c) 
and 3(d). We further calculated the entanglement distribution 
in Eq. (21) and found that the squared negativity is monog¬ 
amous in the composite cavity-reservoir systems. Therefore, 
the quantity can serve as a multipartite entan¬ 

glement indicator as time evolves, as plotted (solid-blue line) 
in Fig. 3 for different probabilities. 

Next, we analyze the MIN distribution of the MMES in the 
multipartite cavity-reservoir system. It has been proved that 
the MIN is not monogamous in multiqubit systems ll4^ Iddll . 
However, whether the MIN is monogamous in multipartite 
multilevel systems needs to be further investigated, especially 
for the newly introduced MMES. Using the relationships of 
the density matrices PciC 2 > Pcir 2 ^ PriC 2 ^ Prir 2 ^ we can get 

MIN,,, 2(0 — [MINe„,(f)], (24) 

where is the exchanging operation acting on the param¬ 
eters ^ and X- After a derivation similar to that for MINciC 2 
we can obtain the MIN of the subsystem Cir 2 , 


-^'^ciri|c2r-2 (0 ~ -^ciri|c2r-2 (^) ^ciC2^^) ^cxr2^'^) 

( 21 ) 

As the system evolves, the four-partite negativity is invari¬ 
ant and satisfies the relation 1 , 3,2 (f) = W,,,, 1 , 3,2 ( 0 ) = 
^, 1 , 2 ( 0 ) = 0.5 since the local dissipation is unitary and the 
two reservoirs are in the vacuum state initially. At a later 
time, the state of the two reservoirs has a form similar to 
that of the two cavities, and we get the relation Prir 2 {t) — 
S^^^[pciC 2 {t)] in which is an operation exchanging 

two parameters {i.e., ^ ^ x x 0. Thus the negativity 
of the reservoirs is 

Nrir2{t) = S'jo-x [.^CiC2 (i)]- (22) 


MIN„, 3 (f) = ^ex^[p^ + 2V3p(l - p)e^ + Fi], (25) 

with the parameter Fi = (30 + 60 x^ + X*)(l ~ P)^- in 
addition, for the MIN of multipartite cavity-reservoir systems 
in the partition Ciri|c 2 r 2 , we can obtain the expression 

MIN„„|, 3 , 3 (f) =MIN„„ 1 , 3 , 3 (0) = i(l-2p+2p2), (26) 

where the MIN is invariant as the time increases because 
the evolution operation (7,,,, {H, t) 0 C/, 2,2 t) is lo¬ 
cally unitary. In addition, we calculate the MIN distri¬ 
bution = MIN„„ 1 , 3,3 (f) - MIN,„ 3 (f) - 

MIN ,,,2 (f) — MIN ,,,2 (t) — MIN ,,,2 (f), which is written as 


We can also derive the relationship p ,,,2 (t) = [Pcir 2 (/)] 
and the negativities 

iV,„ 3 (f) = 54^0iV,„3(f)] (23) 

for subsystems cir 2 and riC 2 . A more detailed description of 
the density matrix Pcir 2 {t) nnd its negativity W,,r 2 (f) can be 
found in Appendix C. 

In Eig. 3, we plot the negativity distribution as a function 
of the time nt for the cases where the probability p of the 
MMES is chosen to be 0,0.5,0.75, and 1. As time increases, 
the two reservoirs exhibit the phenomenon of entanglement 
sudden birth (ESB) ifl^ . which corresponds to the ESD of 
the two cavities. As the probability p increases, the ESB time 
is advanced, and the ESD time is delayed, as shown in Eigs. 
3(a)-3(c). When the probability is p = 1, both the ESB and 
ESD phenomena disappear, as shown in Eig. 3(d), since the 


M^,rPc2r2it) = (1 - V3p), (27) 

with the parameter Gi = (1 — p)(l — X^ + X^)- 

In Eig. 4, we plot the distribution of the MIN as a func¬ 
tion of time Kt with different probabilities p for the MMES. 
As time increases, the MIN of two cavities decreases asymp¬ 
totically, and the MIN of the two reservoirs increases asymp¬ 
totically. When the time nt —?> 00 , the MIN of the two cavi¬ 
ties transfers completely to the reservoirs. Eor the subsystems 
Cir 2 and riC 2 , the MINs first increase to their maximums and 
then decay asymptotically. As the probability increases, the 
distance between the two peaks of MIN ,,,3 and MIN,, ,3 be¬ 
comes smaller. When the probability is p = 1, the distance 
goes to zero, and the two MINs coincide completely, as shown 
in Eig. 4(d). Unlike the distribution of entanglement negativ¬ 
ity, the MIN in the multipartite systems is not always monog¬ 
amous. When the mixed-state probability is p = 0, the MIN 
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FIG. 4: (Color online) The MIN distribution of the MMES in mul¬ 
tipartite 2 (g) 2 (g) 4 (g) 4 cavity-reservoir systems, where the residual 
nonlocality can be positive, zero, or negative as a func¬ 

tion of the time parameter Kt. 


is monogamous and the indicator |c 2 t -2 solid blue 

line) is nonnegative, as shown in Fig. 4(a). However, when the 
probabilities are p = 0.5 and p = 0.75, the MINs are polyg¬ 
amous, and (f) is no longer positive [see Fig. 4(b) 

and 4(c)]. When the probability is p = 1, the MMES becomes 
the two-qubit Bell state, and the indicator (t) is zero 

at all times for cavity-reservoir systems as shown in Fig. 4(d). 
The different distribution property from that of entanglement 
negativity indicates that the MIN and entanglement are two 
inequivalent types of resources for quantum information pro¬ 
cessing. 

Although the MIN itself is not monogamous in multipartite 
multilevel systems, its functions m ay p ossess this property. 
For example, the quantum discord 11^ is not monoga¬ 
mous even in three-qubit pure states 11441 - 1^ . but the squared 
quantum discord is monogamous in an arbitrary three-qubit 
pure state . Recently, similar situations for the entan¬ 

glement of formation have also been discussed ll^ [sol - l^ . 
For the MIN in multipartite cavity-reservoir systems, we cal¬ 
culated the square of the MIN, and the numerical result sup¬ 
ports the monogamy relation. However, in the general case, 
an analytical proof for the monogamy property of the squared 
MIN is still an open problem. 


IV. DISCUSSION AND CONCLUSION 

We have studied the dynamic behavior of the MMES over 
the course of the dissipative evolution of multipartite multi¬ 
level cavity-reservoir systems. It has been found that, unlike 
the two-qubit Bell state \ipi) = (|00)4- IT) )/-\/2, whose neg¬ 
ativity decays in an asymptotic manner Olhll . the entanglement 
dynamics of the MMES exhibits the ESD phenomenon, as 
shown in Eig. 1. Therefore, as an entanglement resource, the 


MMES is not superior to the pure two-qubit Bell state in this 
dissipative system. We think that the high-dimensional com¬ 
ponent 1 ^ 2 ) = (|02) 4- |13))/'\/2 in the MMES gives rise to 
the ESD of the two cavities’ evolution. Moreover, we further 
study the MMESs in 2(g)6 and 2®8 systems where the compo¬ 
nent I' 02 ) is replaced by the higher-dimensional components 
103) = (|04) + |15))/y2 and |04) = (|06) + |17))/V2, re¬ 
spectively. The analytical results show that the new MMESs 
still experience the ESD in the dynamical evolution (the de¬ 
tails for the calculation are given in Appendix D), which fur¬ 
ther supports our viewpoint. 

The MIN has a close relation with quantum communica¬ 
tion protocols involving local measurement and a comparison 
between the pre- and postmeasurement states d. We hnd 
that maximal entanglement cannot guarantee maximal non¬ 
locality. As shown in Eig. 2, the MIN of the MMES is not 
maximal, and its value is directly proportional to the purity 
Tr(p^) of the MMES, which is quite different from the situa¬ 
tion of the Bell state exhibiting maximal nonlocality. Eor the 
MMESs with higher-dimensional components, their MINs are 
also dependent on the mixed state probability p, and the non¬ 
locality evolutions are asymptotical (the details are given in 
Appendix E). We explain that the decrease of the MIN of the 
MMESs is due to the decrease of their purity, and the MIN 
evolution of the MMESs is asymptotic since this kind of non¬ 
locality contains both quantum and classical correlations in. 
Eor the quantum nonlocality related to the violation of Bell in¬ 
equalities ifssl - l^ . its relation to the maximal entangled state 
is still an open problem yet to be addressed. 

In order to obtain a deep understanding of the dynamic 
properties of the MMES, we have investigated its entangle¬ 
ment and nonlocality distributions in multipartite systems. 
The numerical results have shown that the squared negativ¬ 
ity is monogamous in multipartite cavity-reservoir systems 
(beyond the four values of probability p shown in Eig.3, we 
further calculated the distribution forp ranging across [0,1]). 
Moreover, for the MMESs of 2 ® 6 and 2 0 8 systems, our 
numerical calculation still shows that the squared negativity 
is monogamous in the multipartite dissipative systems. These 
results support the conjecture of He and Vidal ll^ that the 
squared negativity is monogamous in multipartite multilevel 
systems. On the other hand, the MIN distribution of the 
MMES is not monogamous in the multipartite cavity-reservoir 
system, as shown in Eig. 4, which indicates that the MIN is a 
different type of resource from entanglement in quantum in¬ 
formation processing. We further investigate the MIN distri¬ 
butions for the MMESs with higher-dimensional components 
and find that the MIN is still not monogamous (the details are 
shown in Appendix E). 

In conclusion, we have studied the dynamic behavior of the 
MMES in multipartite cavity-reservoir systems. It has been 
found that the evolution of the negativity of the MMES ex¬ 
hibits the ESD phenomenon, and is not superior to the two- 
qubit Bell state as an entanglement resource in a dissipative 
system. We also found that maximal entanglement cannot 
guarantee maximal nonlocality. The MIN of the MMES is not 
maximal, and its evolution is dependent on the mixed-state 
probability of the MMES. In addition, we have investigated 
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the distributions of the negativity and the MIN of the MMESs 
in the multipartite cavity-reservoir systems, where two types 
of correlation exhibit different monogamous properties. 
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Appendix A: The derivation of the ESD line for negativity 

In Fq. (12) of the main text, the negativity A"ciC 2 (t) is de¬ 
termined via the sum of absolute values of the negative eigen¬ 
values. After some analysis, we find that the eigenvalues 
{Ai, A 2 , A 4 , Ae, As} are always nonnegative, while the other 
eigenvalues {AsjAsjAvj can be positive, zero, or negative. 
Therefore, as the two cavities evolve, the negativity A"ciC 2 (0 
becomes zero when the three eigenvalues {AsjAsjAzj be¬ 
come nonnegative. 

Using the expressions for A 3 , A 5 , and A 7 in Fq. (13), we 
can derive the p ^ td, relations when these eigenvalues are 
zero. When A 3 = 0, we have the relation 

3(e'"* - 1)2(3-6e'"*-f2e2'‘‘) 

^ ~ 9 - 36e«* -f - 24e3«* -f 2e4«*' 

For the case A 5 = 0, we have 

Atf = ln[(3-f ^)/2] (A2) 

for an arbitrary value of parameter p. When A 7 = 0, we can 
derive the p ^ nt relation as shown in Fq. (14) of the main 
text. In Fig.5, we plot the three relations in the plane of pa¬ 
rameters p and Kt, where the dashed green line is for A 3 = 0, 
the dot-dashed blue line is for A 5 = 0 , and the solid red line 
is for A 7 = 0. The three lines divide the whole area into 
four parts. In regions I, II and III, the signs of the eigenvalues 
(A 3 , As, A 7 ) are (-,-,-), (-f, -, -), and {+, +, -), which 
result in nonzero negativity for the two cavities. In region IV, 
all the eigenvalues are positive leading to the negativity being 
Ncic 2 {t) = 0- Thus, as seen from Fig. 5, the red line for 
A 7 = 0 determines the FSD time of the two cavities, which is 
described by Fq. (14) of the main text. 


Appendix B: Calculation and continuity analysis for the MIN of 
two cavities 

Before evaluating the nonlocality MINcic2(0 given in Fq. 
(18), we first calculate the local Bloch vector x and correlation 



FIG. 5: (Color online) Four regions in the entanglement evolution of 
two cavities, where, in region IV, all three eigenvalues (A 3 , As, A 7 ) 
are positive and the negativity Vcjc 2 becomes zero. 


matrix T. According to Fq. (17), the local Bloch vector of 
subsystem ci is 



which leads to the norm being ||x ||2 = x^/ 8 , with x = 
v/l — The correlation matrix T = T'l2y/2 is a 3 x 15 
matrix, in which the nonzero elements of T' are 

Til = -T2,2=?"[p+(1-p)C" + V3(1-p)x^], 

T',3 = (1-2x^ + 2x")[1-4(1-p)^2^2]^ 

TU = Tiiii = nfi = -Tl^ = VQ{l-p)ex\ 

Ti',i3 = + 

Tin = -^^[P - - P)i^ + Vsil - p)x% 

T 3.12 = X^-4(l-p)eV, 

r'_i 5 = 2p - 1 + ( 6 p - 4)x" - (8 - 10p)x" 

+6(1-p)x®, (B2) 

with ^ When nt > 0, we have the local Bloch 

vectorx 7 ^ 0. After substituting the three terms x, | |x| p, and T 
into the first formula in Fq. (18), we can obtain the expression 
forMINciC 2 (^ > 0) in Fq. (19). When Kt = 0, the quantum 
state /OciC 2 ( 0 ) is the MMES for which the local Bloch vector 
is X = 0. In this case, we need to calculate the eigenvalues of 
matrix TT*, which are Ai = A 2 = A 3 = (1 — 2p 4 - 2p2)/4. 
According to the second formula in Fq. (18), we can derive 
MIN,,,,(0) = (l-2p + 2p2)/2. 

Next, we prove the continuity of MINciC 2 (t)- Based on the 
previous analysis, we know that the two formulas in Fq. (18) 
are used for the cases x 7 ^ 0 and x = 0 , which correspond to 
the time evolutions Kt > 0 and Kt = 0, respectively. If the 
MINciC 2 is continuous, the limit of MINciC 2 (ttf —f 0 +) in the 
first formula should coincide with the value of MINciC 2 (ttf = 









0). After some calculation, we can get 


lim 

Kt—¥ 0 ^ 


x*Tr*x 


lim 

K.t^ 0 + 


d(x*TT‘x) 

d{Kt) 

4l|xP) 

d{Kt) 


lim 

Kt—¥ 0 + 


d 2 

d{Kt)‘^ 

d(Kty^ 


±{l- 2 p + 2 p^) 

1 

4 

^3, 


(B3) 


where we have used L’Hopital’s rule and A 3 is the minimal 
eigenvalue of matrix TT*. Then, in the limit nt 0+, the 
two formulas in Eq. (18) are continuous and we have 

lim = MINcic 2 ( 0 ) 

= i(l- 2 p + 2p2). (B4) 

As a result, the nonlocality MINciC 2 (^) can be described by 
Eq. (19) of the main text throughout the entire period of the 
dynamic evolution. 


Appendix C: The density matrix Pcir2 and its negativity A^cir2 


Throughout the dynamic evolution of multipartite cavity- 
reservoir systems, the quantum state of subsystem cir2 is 
Pcir2 (t) = trr-iC2 [PciriC2r2 (^)]> which Can be written as 


Pcir2 (t) 


/ 611 0 0 0 0 

0 622 0 0 0 

0 0 633 0 0 

0 0 0 644 0 

0 0 0 0 655 

661 0 0 0 0 

0 672 0 0 0 

Vo 0 633 0 0 


6i6 0 0 \ 

0 627 0 

0 0 638 

0 0 0 
0 0 0 
bee 0 0 

0 677 0 

0 0 bs8 / 


, (Cl) 


where the nonzero matrix elements are 


6 n = [p+(l-p)a(l + CV)A 

622 = {2(1 - + [p + 3(1 - p) Ax'‘}/2, 

bes = (1 - P)x‘‘(l + 3CA^)/2, 

644 = (1 -p)x®/ 2 , 

A = F+p(e4-A]/2, 

bee = + 3(1 - p) A/2, 

677 = 3(1 -p)^‘‘x'‘/2, 

688 = (1 -p)?V/2, 

6i6 = bei = ^x[p + v^(l - ?')A/2, 

627 = bf 2 = \/3/2(l — p)C^x^, 

&38 = &83 = (1-p)AV2- (C2) 

Eor this quantum state, the negativity is 


8 

where Ai are the eigenvalues of the partial transpose matrix 

T 

Pcih and have the form 

Ai = (1-pAxV2, 

A2 = (i + ?V)b + (i-p)A/2, 

A 3 = [1 - 2 x" + Bi - A-(1 -p)x"S2]/4, 

A4 = [1 - 2 x" + Si + V1-(1 -p)x2 S2]/4, 

As = [(3 - 2p)x^ + Be- y/F^]/4, 

Ae = [(3 - 2p)x" + S 3 + \/F^]/4, 

At = [(1 - p)x'‘(3C‘‘ + X*) - \/(l -p)2S5]/4, 

As = [(1 -p)x"(3eAx") + V(1 -p)"A]/4,(C4) 

with the parameters 

Si = (7 - 5p)x^ - 10(1 - p)x® + (4 - 4p)x®, 

B 2 = (16 - 8V3)p + [14 - 24(2 - V3)p]x^ 

- 8[8 - (10 - 3V3)p]x^ + [123 - (111 - 8V3)p]x® 
-(104 - 96p)x® + 28(1 - p)x'° + (8 - 8 p)x'" 
-4(1-p)x^", 

5 3 = -( 8 -7p)x‘‘+ 12 ( 1 -p)x® - 6 ( 1 -p)x®, 

5 4 = (3 - 2pf - [36 - 2(23 - 6 p)p]x^ 

+(64 — 96p + 33p^)x^ — 12(1 — p)(4 — 3p)x® 
-12(1 -p)A®, 

Ss = X®(9f+ 4eV-6eV + X®)- (C5) 


Appendix D: The ESD for the MMESs with higher-dimensional 
components 

In the multipartite cavity-reservoir systems, we first con¬ 
sider that the two cavities are initially in the MMES 

AL(0) =p|bl)(bl| + (1 -p)lb3)(b3|, (Dl) 

where l-i/ji) = (|00) + \ll))/y/2 is the two-qubit Bell state 
and \'ij^e) = (|04) + |15))/-\/2 is the high dimensional com¬ 
ponent. Along with the evolution of cavity-reservoir systems, 
the output state is 

PciVc27'2(^) — ly\i\'t’o)c^r^\4’o)c2r2 ~‘<~\4'l)c^r^\4’l)c2r2) 

Am ciri (M C2r2 C2l^2)] 

1 — f) 

+ ciri\4'4) C2r2 4“ bl) ciri j </5)c2r2 ) 

x((</o| ciri (04|c2r2 + (</l|ciri(</5l C2r2 )]CD2) 

where the components have the forms 

Ao) = |00), 

\A) = /|io) + x|oi), 

lA) = /"|40) + 2^A|31) + V6eV|22) 

+2exA3) + xl04), 

14) = f|50) + 45eA|41) + 4T04x"|32) 

4l0CV|23) +45^X^114)+x"|05), (D3) 


lVcir-2 (t) 


2 



9 




has the same form as that in Eq. (D2), but the components 
and \(f>\) are replaced by the new components \(j)\) and 
\4 >t), which can be written as 

14) = e®|60) + 46e"x|51) + 4T5CV|42) 

+4^4 + V6ex^|15) 

+X®|06), 

14 ) = 4 | 70 ) + 47 e®x| 61 ) + 4 ^CV| 52 ) 

+\/^4x^| 43) + 4^4x134) + 4^4x^|25) 

+475x®|16) + x"|07). (D7) 


FIG. 6: (Color online) The BSD lines in the evolution of two cavities 
which are initially in the MMESs: (a) ™ Eq. (Dl) and (b) piic2 

in Eq. (D6). 


with the parameters ^(f) = e and x(f) = (1 — e 
By tracing the subsystems of two reservoirs, we can get the 

output state of two cavities PciC 2 (i) = Trrir 2 [PciriC 2 r 2 (f)]. 
which is a 12 X 12 matrix. In order to obtain the entangle¬ 
ment negativity of we calculate the eigenvalues of 

(l)Tc 

the partial transpose matrix Pci^ {t). After some derivation, 
we find that there are four eigenvalues which can be negative. 


A 2 = (l-p) 4 °( 3 x"-v^r+ 44 )/ 2 , 

Ae = (1 - + 15x^ - 4l + 70x4 + 25x8)/4, 

At = (l-p)eV(l + 5x"-41 + 15x4, 

All = 4x4x42 +3x4+P(l-2x"-3x®)]/2 

-4^/2, (D4) 

where the parameter is iJi = + 245(1 — p)px^ + 

(1 — p)^(4x^^ + 13x^® + 4x^°)]. Similar to the analysis 
in Appendix A, we can derive the ESD time for the MMES 
according to the four eigenvalues. When the mixed- 
state probability p changes in the region [0,pi] with pi = 
(347 - 12545)/1922 r; 0.03512, the ESD time for the 
MMES is Kt = ln[(5 -F 45)/4] « 0.5928. When the proba¬ 
bility p G [pi, 1), the ESD time is determined by the p ^ Kt 
relation 


_ 42 x ^44 + x "+2 _ 

1 — X"^ + 2(2 — 45)x® + 6x^^ + X^® ~ 5x^'’ 


(D5) 


where the two parameters are Ji = 4 — 245 -F 3(3 — 45)x‘^ + 
2x® and J 2 = 2 — 45 -F 3x"^ + X* ~ 5x^^. In Fig.6(a), 
we plot the ESD line (red line) as a function p{Kt), which 
divides the entanglement evolution into an entangled region 
and a disentangled region. 

Next, we consider the two cavities which are initially in the 
MMES, 

^444) =plV'i)(4l + (1 -p)|V'4)(V’4|, (D6) 


where |■^/> 4 ) = (|06) -F 117))/42 is the high-dimensional com- 

(o') 

ponent. As the systems evolves, the output state PciriC 2 r 2 (1) 


After tracing the subsystems rir 2 , we can get the output state 
f 2) 

of two cavities pci^t). Furthermore, by doing the partial 

(2)Tc 

transposition, we can obtain the matrix pciC 2 ^ {t) and calcu¬ 
late its eigenvalues. The ESD line is determined by the neg- 

(2)T 

ative eigenvalues of Pci^ (t). When the mixed state prob¬ 
ability p G [0,p2], with p 2 = (8669 — 240l47)/370191 « 
0.006258, the ESD occurs at the time Kt = ln[(7 -F 47)/6] ~ 
0.4748. When p G [p 2 ,1), the ESD time is determined by the 
following p ^ Kt relation: 


_ x^HKi + yi^) _ 

1 - x^ + 2(3 - 47 )x^ 2 + 8x^6 + X^^ - 7x2® 


(D8) 


where the two parameters are ATi = 3—47+4x^+X^^~7x^® 
and iXa = 15 - 647 -F 8(4 - 47)x'^ + 9x®. In Fig.6(b), 
we plot the ESD line (blue line) as a function p{Kt), which 
cut the entanglement evolution region into two parts, i.e., an 
entangled region and a disentangled region. 


Appendix E: The MIN of the MMES with higher-dimensional 
components and its distribution 

We first consider the MMES piiC2(0) with the high¬ 
dimensional component I'i/a) = (|04) -F |15))/42, as shown 
in Eq. (Dl). According to the formula in Eq. (18) of the main 
text, we can derive the nonlocality 

(0)] = {p- 1/2)2 ^ (Ej) 

which is dependent on the mixed-state probability p and di¬ 
rectly proportional to the purity of the MMES. In the cal¬ 
culation of the MIN, the matrix basis for the subsystem C 2 
is chosen to be the generalized Gell-Mann matrices (GGM) 
which are the higher-dimensional extension of the Pauli 
matrices. The GGM basis for a d-dimensional system is com¬ 
posed of three types of matrices IHl]: (i) d{d—l)/2 symmetric 
GGM, 


Af = |j)(fc| + |fc)(j|, l<j<fc<d; (E2) 
(ii) d{d — l)/2 antisymmetric GGM, 

Af = -Aj){k\ + i\k){j\, l<j <k<d; (E3) 
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FIG. 7: (Color online) The MIN distributions of the MMESs in mul¬ 
tipartite 2 ® 2 (g) 6 (8) 6 and 2 (g) 2 (g) 8 (g) 8 cavity-reservoir systems are 
plotted as a function of Kt, where the negative values indicate that 
the MIN is not monogamous. 


and (iii) (d — 1) diagonal GGM, 

i 

A'= v'2/((2 - |j)(j|_/|) + l)(/ + l|) (E4) 

j=i 

with 1 < Z < d — 1. It should be noted that the GGM needs to 
be normalized in the generalized Bloch form of /oiiC 2 • Along 
with the interaction between the cavities and reservoirs, the 
MMES will evolve into a 2 0 6 system. After some derivation, 
we obtain 

MIN[p«,^(f)] = ^l{£}^+p[p-{2-p)e^]+L^+L2} (E5) 

where the two parameters are Li = 2(1 — p)[^/bp + 30(1 — 
and L 2 = (1 — In 

the dissipative procedure of cavity-reservoir systems, the non¬ 
locality of two cavities decays in an asymptotical way, which 
is similar to the case of MMES in 2 (g) 4 systems. 


Next, we analyze the MMES of a 2 (g) 8 system in Eq. (D 6 ), 
which has the high-dimensional component = (|06) + 
|17))/V2. It is found that the MIN for this MMES pi?c 2 ( 0 ) 
has the same expression as that in Eq. (El), which is also 
dependent on the mixed-state probability p. As the system 
evolves, the MIN for two cavities decays asymptotically and 
can be expressed as 

MIN[pg) Jf)] = ^[^24 + L 3 + {l-pfU], (E 6 ) 

where the parameters are L 3 = p[p — (2 — + 2(1 — 

p)[-\/7p+350(l—andL 4 = 42^^°x‘^ + 315^^®x® + 

525C®x^® + 126C'‘x^° + 

Eor the MMESs piiC 2 ( 0 ) and with the high¬ 

dimensional components, we further calculate the distribution 

|c 2 r 2 (^) “ MIN(,jri|c 2 r -2 (0 ~ ^II'IciC 2 (i) 
-MIN,,,,(f)-MIN,,,,(f) 
-MIN,,,,(f) (E7) 

in the multipartite cavity-reservoir systems. We find that the 
MIN distributions are still not monogamous. As examples, 
we choose the mixed-state probability p = 0.8 for the two 
MMESs and calculate their MIN distributions. In Eig. 7, 
the MIN distributions in the multipartite systems are plotted 
, where two cavities are initially in the MMESs (0) and 

PciC 2 ( 0 ). As shown, the negative values for the distributions 
indicate that the MIN is not monogamous. 

However, for the squared negativity of the MMESs in 2 ® 6 
and 2 g) 8 systems, we calculate the entanglement distribution 
in the multipartite 2 g 2 g 6 g 6 and 2 g 2 g 8 g 8 cavity- 
reservoir systems, where the mixed-state probability p ranges 
across [0,1]. The numerical results still support that the nega¬ 
tivity is monogamous. 
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